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The energy for the coronal heating must be provided from the convection zone. 
However, the amount and the method by which this energy is transferred into the 
corona depends on the properties of the lower atmosphere and the corona itself. 
We review: 1) how the energy could be built in the lower solar atmosphere; 2) how 
this energy is transferred through the solar atmosphere; and 3) how the energy is 
finally dissipated in the chromosphere and/or corona. Any mechanism of energy 
transport has to deal with the various physical processes in the lower atmosphere. 
We will focus on a physical process that seems to be highly important in the chro¬ 
mosphere and not deeply studied until recently: the ion-neutral interaction effects 
in the chromosphere. We review the relevance and the role of the partial ionization 
in the chromosphere and show that this process actually impacts considerably the 
outer solar atmospher e. We include analysis of our 2.5D radiative MHD simulations 
with the Bifrost code ( Gudiksen et al. 201lll including the partial ionization effects 
on the chromosphere and corona and thermal conduction along magnetic field lines. 
The photosphere, chromosphere and transition region are partially ionized and the 
interaction between ionized particles and neutral particles has important conse¬ 
quences on the magneto-thermodynamics of these layers. The partial ionization 
effects are treated using generalized Ohm’s law, i.e., we consider the Hall term and 
the ambipolar diffusion (Pedersen dissipation) in the induction equation. The in¬ 
teraction between the different species affects the modeled atmosphere as follows: 
1) the ambipolar diffusion dissipates magnetic energy and increases the minimum 
temperature in the chromosphere; 2) the upper chromosphere may get heated and 
expanded over a greater range of heights. These processes reveal appreciable differ¬ 
ences between the modeled atmospheres of simulations with and without ion-neutral 
interaction effects. 

Keywords: Magnetohydrodynamics MHD —Methods: numerical — Sun: 
atmosphere — Sun: magnetic field 


1. Coronal heating vs heating concentrated in the low 

atmosphere 

One of the recent and heated discussions on chromospheric and corona heating 
concerns the location of magnetic energy release in the solar atmosphere. Is it 
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highly dependent on height? What is the role of flux emergence and what is the 
role of the magnetic field topology? 

One dim ensional hyd r odyna mic simulations can cast some light on this issue. 
For example, Testa etal\ ( 201411 use ID models to prove that nano-flares situated in 
the corona can spread energy via non-the rmal electron beams rather th an through 
thermal conduction. As another example, Klimchuk fc Bradshaw ( 2014[) compared 
synthetic data obtained from ID multi parametr ic hydro simulations with the At¬ 
mosp heric Imaging Assembly observations (AIA, iLemen et al [Boeri^r et a,I 


1201^ on board of the Solar Dynamics Observatory ('SDO. IPesnell et a/.ll2012ll . They 
claimed that their synthetic data can not reproduce the observables unless they 
consider quasi-steady coronal heating or coronal nano flares situated well above the 
chromosphere. This particular case highlights the limits of a ID hydrodynamic ap¬ 
proach: advanced 3D simulations of heating in association with chromospheric jets 
indicate that the complex 3D topology needs to b e taken into account to understand 
the impact on the corona ( Hansteen et a,I ][2^. 

It is clear that studying the emergence and topology of the magnetic field and the 
likely sites of energy release requires full 2D or 3D models covering the region from 
the photosphere to the corona. So phisticated simulations that in clude thermal con¬ 
duction along magnetic field li nes ( Gudiksen fc Nordhmd 2004h . radiative losses in 
the corona and chromosphere ( Hansteen et Qi.l[2007 1 are now well e stablished tools 


and provide quite sophistic ated comparisons with observations (IHansteen et al 


and provide quite sop 
20101 : Bingert fc Peteii 1 


20111) . These models show that the average heating per par¬ 


ticle coming from the released magnetic energy is largely concentrated in the region 
between the upper chromosphere and the lower corona. Synthetic observations gen- 


shifts, widths, and profiles (see e.g.. 

Gudiksen & Nordhmdl 20056: Peter et oi.l2006: 

Olluri et al. 201361: 

Hansteen et all 

2014^. For instance. Hansteen et al. (20101) use 

the Bifrost code (Gudiksen et al. 2011) and achieve rather good agreement between 


synthetic and observed average redshifts for transition region spectral lines. They 
argue that this is due to episodic heating being localized in the lower atmosphere 
(somewhere betwe en the middle chromosphere and the lower corona). 

Simulations by Martfnez-Svkora et all ( 2011al) also reveal another observable 


due to the presence of localized episodic heating: they studied the asymmetries 
of the EUV spectral pr ofiles using so-called Red-Blue asymmetry (RB) analysis 


gion EUV profiles (Kieldseth Moe & NicolasI 

1977: Peter 

2001 

Hara et al 

200(4 

De Pontien et aZ.ll2009 

: McIntosh & De PontieiJ2009l: IPetei 

2010 

:lTian et a,l. 

2011). 


lines from their simulations reproduces, qualitatively, the obse rvations taken by So- 


lar Ultraviolet Measurements of Emitted Radiation (SUMER JWilheln^^t_al 


on board of the Solar and Heliospheric O bservatory (SOHO. IPomingo et al. 


199.U 




and the EUV Imagi ng Spectrometer ('ElS. ICnlhane et fli.ll200(Tl) on board of Hinode 


(IKosugi et a/.ll200lt) . Eurthermore, they find that the episodic heating produces a 
larger range of velo cities within the transit i on regio n than in the corona. Moreover, 
the simulations by Martfnez-Svkora et all ( 2011a ) include small scale flux emer¬ 
gence, and find that t he average volumetric heating is still localized large ly within 
the transition region ( Martfnez-Svkora et o,l\\20()9a : Hansteen et ai.lEoiOl) . Despite 
this, they find a quantitative disagreement with observations when looking at line 
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asymmetries, most likely due to the lack of vigorous dynamics in the simulations. 
A word of caution is in order, these synthetic spectral results and models may be 
significantly impacted by the simplified magnetic field configuration utilized in the 
simulations. The penalty paid in performing full 3D simulations is the high resis¬ 
tivity required to keep structures larger than the numerical grid scale. Small scale 
plasma physics can therefore not be fully resolved, and the treatment of phenom¬ 
ena such as the damping of shocks and the magnetic reconnection process occur at 
scales much larger than in nature. Therefore, one must be aware of the limitations 
of not resolving some of the physical processes which may impact the transport of 
the magnetic field stresses (built by the convective motions) to greater heights in 
the upper atmosphere. 


2. How to build stress in the atmosphere 


Let us start at the source: it is well known and broadly accepted that convec¬ 
tive motions build stresses in the magnetic field in the atmosphere. In the quiet sun 
and coronal holes, there are at least five different mechanisms in the photosphere 
that could build stress or transport energy into the upper layers. The first is the 


ing. convective collapse, etc (Carlsson & Stein! 1992. 1995: Bellot Rubio et oi.!2001: 

Stein & Nordlundl 2006: Skartlien et all 

2000l:!Hanstee 

n et all 2006!: De Pontieu et al. 

2007a!: Martinez-Svkora et a7!!20096l: Kato et a/.!!2011 

among many others), or even 


by reconnection or magnetic energy release in the photosphere or lower chromo¬ 


sphere processes (jMartfnez-Svkora et aLN20096ll . The second lies in waves, where 
the restoring force is “magnetic” such as kink-, Aflven-, or other transversal waves. 


vective motions (De Pontieu et all 200761: !.less et al 20091: iGoossens et all 201^ 

McIntosh et aLll2011: van Balleffooiien et a/Jl201ll: De Moortel & Pascoel2012 

!De Pontieu et all 

2014 

). The third is by small scale flux emergence (Centeno et al. 

2 OO 7 I: Martinez Gonzalez & Bellot Rubiol 

2009 

: Martinez Gonzalez et al. 2012: Sainz Dalda et al. 2012 

among others). In 


short, these ob servations reveal that sm a ll sca le loops occur throughout the pho¬ 
tosphere, e.g., iMartfnez Gonzalez et all (1201211 observations of the quiet sun re¬ 


vealed the emergence of small scale magnetic field loops in the photosphere at 


a rate of the order of 0.2 loops h“^ arcsec“^. The fourth mechanism is through 
mass flux. The proponents note that mass flux also can contribute to transport 
or build energy in the upp er layers of the atm osphere, for instance photospheric 
supers onic jets observed by Bellot Rubiol (|2009ll with t he Solar Optic a l Tele scope 
fSOT. IdAuneta ef^. 2008) on board of Hinode and bv iBorrero et al\ (2012) with 
IMAX on board of the balloon-borne solar observatory SUNRISE ( Solanki et all 
I2OIOI) . Another example of mass flux contributing to energy supply in the outer 
layers is that of spicules, with De Pontieu et all ( 201lll suggesting that some of the 
spicular mass flux is heated to coron al temperatures, based on observat ions with 
Hinode/SOT and SDO/AIA (see also Martfnez-Svkora et a,l. 2011611201^ . Finally, 
the motion of magnetic flux elements in the photosphere, driven by convective mo¬ 
tions will easily build magnetic field stres ses tha t prop agate into the corona, as 
illustrated with the well known cartoon by Parker! ( 1983ll . Note, these mechanisms 
are not completely independent, e.g., flux emergence can also drive waves. Another 
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possible source of the mechanical energy needed to build extra current is the highly 
dynamic chromosphere (see below). 

We would like to describe in greater detail two of the sources listed above that 
we consider of great interest in this review. The first one concerns small scale flux 
emergence. Small bipole structures in t he Stokes V magne t ograms suggest tin^y loops 

within the g ranules such as obser v ed by Lites et ^ 11996 ^ 1998ll. De Pontieu J 2002 j _ 

and la ter by Centeno_^t_a^ (2007), Ishikawa et al. ~()2008 UMartmez Gonzalez fc Bellot Rubio 
( 2009ll . and iGbmorv et aJ ( 20inll . among others. A large number of small scale loops 
are found to occur in the quiet sun. In addition, stro ng asymmetric Stokes V profiles 
suggest structures of unresolved, even t i nier, loons (iGrossruann-Doerth et al. 2000l : 

Sainz Dalda et al. 201‘A Viticchie 2012 1. Sainz Dalda et al. ( 2012 1 combined inver¬ 


sion with synthetic profiles fro m 3D radiative MHD s imulations in order to explain 
the strong asymmetry profiles. Viticchie et all ( 201lll categorized a large sample of 


different asymmetries and inverted the corresponding Stokes profiles. Many of these 
strong asymmetries suggest the presence of two different polarities of magnetic field 


amount of horizontal magnetic flux at the ohotosohere (ILites et al. 

19961 2008: 

iBorrero & Kobell 2011: 

Bellot Rubio & Orozco Suarez 2012; Stenflo 

2 OI 3 I among 

othersl. IMartfnez Gonzalez et al. 

(2012 

) suggested that 30% of these tiny loops are 


able to reach the chromosphere (observa t ions, for instance, by Centeno.^t_aLl 


iMartfnez Gonzalez fc Bellot Rubi3 20091 : iGuglielmino et '5]l201oi : l0rtiz et al. 


2007 


120141 


also reveal the various chromospheric signals of flux emergence). 

The second aspect we would like to discuss in a bit more detail is the mo¬ 
tion of the magnetic flux elements. Magnetic field elements move around due to 
the convective motions, collapsing granules and emerging new convective cells. The 
magnetic field elements evolve with the convective motions, merging elements of the 
same polarity, fra gmenting magnetic e lements, and canceling magnetic elements of 
different polarity (ISchriiver et aLlll997l) . It is important to understand how the con- 
vective motion is braiding and tangling the magnetic held in the solar atmosphere 
( Lamb et al.ll200^ IPeForest et al. 20071 : Zhou et al. 2010l Giannattasio et al. 20141 
among others) since the convective motion could build magnetic energy and there¬ 
after this energy may dissipate into kinetic and thermal energy. These authors 
studied the dynamic properties of the magnetic field elements at the photosphere 
in the quiet sun. Most of the identified magnetic network elements seem to form due 
to unresolved small scale shredding, coalescence of previous magnetic flux, magnetic 
field cancelation and emergence. 

Moreover, it is crucial to calculate the magnetic elements that have merged, can¬ 
celed and/or fra gmented in order to understand how the magnetic field i s stressed in 


celed and/or fra gmented m order to understand now tde magnetic neld i s stressed ir 
the a t mosphere (ISchriiver et "^119971 llida et a]]l2012t IZhou et ^l2013l iGosic et~d. 


2 OI 4 II . iGosic et al\ (I 2 OI 4 I ) using long time series of Hinode/SOT observations fol¬ 
lowed and identified different internetwork elements, and noting which of these 
reach the network. Their conclusions are that internetwork elements interact with 
network elements and thereby modify its net magnetic flux. Only a tiny part of the 
internetwork elements that interact with the network comes from flux emergence 
and most of the network flux comes from shredding and coalescence of previous 
internetwork magnetic flux. The internetwork magnetic flux replaces the network 
flux in a time period of 9 hours. On average . 40% of the inter network individual 
elements ends up, eventually, in the network ( Gosic et al. 2014[l . 
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3. The various “filters” of the various atmospheric layers 


As discussed above, convective motions have several methods of generating 
enough energy to heat the corona. The next question becomes: What happens 
to this energy for each source mentioned in the previous section as it propagates 
through the solar atmosphere? In the solar atmosphere there are many physical 
processes that play different roles in different layers. Likewise, one must also con¬ 
sider whether the transport of energy can trigger other mechanisms that contribute 
further to increasing the magnetic stresses in the atmosphere. Finally, there remains 
the question of where magnetic energy is released for each source, if it is released 
at all. 

In short, in order to reach the corona, the energy must go through: 

• The photosphere, where a large amount of the convective energy is released 
by radiative transport, and the atmosphere is highly stratified with a pres¬ 
sure scale height of the order of less than 100 km. The upper photosphere is 
dominated by large velocities due to convective overshoot. 


• The chromosphere is the interface layer between the solar surface (photo¬ 
sphere) and the million degree corona. The chromosphere is of great interest 
because, in principle, it contains enough non-thermal energy to heat the en¬ 
tire corona. Many complex physical processes play a role in the chromosphere: 
a) the plasma is in non-local thermodynamic equilibrium (NLTE), b) the ra¬ 
diation is optically thick, c) radiation is highly scattered, d) ionization is 
not necessarily in equilibrium, e) the gas is partially ionized, f) in the upper 
chromosphere, transition region and corona, the thermal conductive flux prop¬ 
agates along the magnetic field lines. In addition, the chromosphere is also 
distinguished by several transitions such as from non-magnetized to magne¬ 
tized, and from partial to full ionization, from high to low plasma /3 (where 
plasma /3 = is the ratio between the gas pressure and magnetic pressure), 
from optically thick to optically thin, etc. Note that many of the processes 
(especially from the items a to d) lead to highly complex interpret ations of 
imaging and spectral observations le.g.. iLeenaarts et a/.ll2012U20I,‘l^ . 


• The transition region, located between the chromosphere and corona, is op¬ 
tically thin, magnetically dominated, with steep temperature gradients, filled 
with shocks and flows approaching or exceeding the local speed of sound, and 
where, as mentioned, thermal conduction is a highly important process. In 
this region we cannot expect ions to be in ionization equilibrium and for low 
density regions such as coronal holes we may also find that ion and electron 
temperatures differ ( Hansteen et al. I993t Lie-Svendsen et all 2001 1. 


• Finally, in the corona, which generally has low plasma beta, radiation is opti¬ 
cally thin, and, as i n the transition region, the ionization of some heavy ions is 


caily thin, and, as m the transition region, the ionization ol some heavy ions is 
not in equilibrium djoselvn et a/.l[l977[ iHansteen et aLl[l9^ iBradshaw &: Car^ 


2006t Olluri et aLll2013a l. background radiation heats lower layers, the timescales 

of several important processes are very short and non-thermal physics may 
be vital in describing these processes. 
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Any kind of magnetic stress built by the convective motions that reaches the 
corona will be strongly affected by the various processes mentioned above. For 
instance, the chromosphere appears to be filled with magnetic field, even in the 
quiet sun, despite the fact that the photosphere is sub-adiabatic, which thwarts 
the expansion of em erging magnetic flux into the regions above ( Acheson 19791 : 
Archontis et al.lboOd l. Therefore, it remains unknown under which conditions small- 


scale flux emerges into the chromosphere. The challenge facing the researcher is thus 
to combine all these ingredients into a model of the entire solar atmosphere, taking 
into account all important effects while discarding those that have little impact on 
the buildup, transport and dissipation of the energy flux heating the outer atmo¬ 
sphere. 


4. Partial ionization effects 


In this paper we will concentrate on ion-neutral interaction effects among all the 
possible physical processes from the photosphere to the lower corona. The chromo¬ 
sphere is partially ionized: ions are coupled to the magnetic field, whereas neutrals 
are not directly affected by the presence of the magnetic field and they can move 
“freely”. Sufficient collisions between ions and neutrals will couple the neutrals to 
the magnetic field, while at the same time to some extent allowing ions to slip 
across the field. As a result, the magnetic field can diffuse and magnetic energy will 
be dissipated into thermal energy. In order to model this scenario, the fluid should 
be treated as consisting of two or three fluids, i.e., ion, neutral and/or electron 
particles. To solve the equations for all three species is computationally highly de- 
manding. Therefore two or three fluid codes focus on spec i fic and localized pr oblems 
( Smith fc Sakail 120081 : LSakai fc Smith 20091: Meier 2011 ; Leake et all 20121 among 
others) instead of modeling large portions of the solar atmosphere, i.e., instead of 
including the convection zone, photosphere, chromosphere, transition region and 
corona. Additionally these codes rarely consider many of the other important phys¬ 
ical processes in the solar atmosphere such as radiative transfer, thermal conduction 
along the magnetic field lines, etc, which is the matter of interest we wish to cover 
in this review (see previous sections). 


(a) The Hall Term and Ambipolar Diffusion 


Fortunately, the chromosphere seems to be sufficiently collisional to allow us 
to convert the three fluid problem into a single MHD fluid while still taking into 
account the most important partial ionization effects. In order to do this one as¬ 
sumes that 1) collisions are sufficiently numerous and the ion and electron temper¬ 
atures are the same, 2) the relevant resolvable timescales one wan ts are assumed 
to be greater than the collision times between ions and neutrals dCowling 1957 : 
Braginskii 19651: Parkeij[2007 : Pandev fc Wardle 2008 : Leake et al. 2013 oil . Accord¬ 


ing to the comparison of the collision rates with timescales of various relevant phys- 
ica l processes from the sel f -consi stent 2D radiative-MHD atmospheric models done 
bv lMartfnez-Svkora et all ( 2012h using the Bifrost code ( Gudiksen et aLll201lll . the 
approximations mentioned above seem to be fulfilled in the chromosphere most of 
the time. 
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In this case, the induction equation expands into the so-called generalized in¬ 
duction equation which includes at least two new ter ms, i.e., the Ha l l term and 
the a mbipolar diffusion (or Pedersen dissipation) term ( Cowlinel Il957l : iBraginskil 
1965) as follows: 


9B 


— = V X 


u X B - r^ohmS - xB + ^(JxB)xB 


\B\ 


(4.1) 


where B, J, u, and rjohm are magnetic field, current density, velocity field, and the 
ohmic diffusion, respectively. The new “diffusion” terms are: 


Vohm — 

'Hhall — 

Vamb — 


meiVei + Ven) 


line 

\B\ 

Qene 

mpjpf {\B\pjpr 


Pi^ii 


Pn^n 


(4.2) 

(4.3) 

(4.4) 


where Pi, Pn, P, t^ni, ne, and are ion density, neutral density, total density, ion- 
neutral collision frequency, neutral-ion collision frequency, electron number density 
and electron charge, respectively. Note that rjhaii is not a diffusive or dissipative term 
in contrast to rjohm and rjamb- The ambipolar diffusion should not be confused with 
the so-called ambipolar drift used in plasma physics. For plasma physics, ambipolar 
drift or electric field is du e to the relative veloc it y or t emperature difference between 
ions and electrons, e.g., Arefiev fc BreizmanI (|2008h . Ambipolar diffusion is also 
referred to as Pedersen dissipation. 

In short, the ohmic diffusion can change directly the magnetic field topology, 
i.e., change the connectivity of the field lines with reconnection, but the Hall and 
the ambipolar terms cannot, without the help of ohmic diffusi on, since the former 
term does not follow the formalism of u x B (iBiskamnI l2005h and the other two 
follow it, where u for the Hall term is J and for the ambipolar term is J x B. The 
ohmic (first term in the right hand side of the following expression) and ambipo¬ 
lar diffusion (second term in the right hand side of the following expression) can 
dissipate magnetic energy into thermal energy as follows: 


de . |J X Bp 

^ Vohml 4 “ Pamb 1^12 


(4.5) 


Isee |Park^l2007l for details). 

The physical process behind the Hall term is that in a fully ionized plasma as in 
the corona, the relative immobility of ions relative to electrons in response to the 
Lorentz force leads to an electric field parallel to J x B. The Hall term for weakly 
ionized plasma is created due to collisions between ions and neutrals. The neutrals 
decouple ions from the magnetic field lines leading again to an electric field parallel 
to J X B. The physical process behind the ambipolar diffusion is that neutrals do 
not experience the Lorentz force. As a result, neutrals decouple from th e magnetic 
field and allow the magnetic field to diffuse through the neutral gas (IBraginskii 
196,5ll . 
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(6) Impact of the Ion-Neutral Interaction Effects on Simplified Physical Processes 


Partial ionization effects impact various physical processes that play a role in the 
chromosphere and prominences. For instance, the partial ionization effects in the 
chromosphere are known to allow ion-neutral col lisions to dissipate alfvenic waves 
( De Pontieu fc Haerendell 1998 : De Pontieu 1999ll . However, the picture is far from 
complete, sinc e it is unknown i f this dissipation plays a major role in heating the 
chromosphere ( BiermanrJ 19^ ^ heak^etji^ 2005 : Hasan fc van Ballegooiien 2008 


Song fc Vasvliilnasll2011 i TMadsen et aLl2014 ). Most of these studies obtain different 

high frequency cut-offs for alfvenic waves for different features (spicules, filaments, 
prominences, thin flux tubes, etc) using fixed values for the io n-neutral collision 


frequency of semi-empirical models such as the VAL-C model (|De Pontieu et al 


2001 : Leake et all 2005t Forteza et all 2007 1. This dissipation contributes to the 


thermal and dynamic energy budget. In fact, as mentioned in the previous section, 
electrical currents perpendicular to the magnetic field can be dis sipated by ambipo- 
lar diffusion and convert magnetic energy into thermal energy ( ArbM_ef_^ 2009 ■ 
Goodman 2011 : Khomenko fc Colladosll201^ Goodman fc Judge l2012ll . Soler et al 


( 20091 . 12 OI 2 I 1 deduced that kink waves are damped by diffusion parallel to the mag¬ 


netic field and Alfven waves by perpendicular diffusion. For short wavelengths, kink 
waves are converted into Alfven waves due to the collisions between ion and neu¬ 
trals. Ohmic dissipation dominates for long wavelengths. In fact, all the diffusive 
mechanisms must be considered , and missing any of the inertial, Hall and ambipolar 
term s will give incorrect results ( Khodachenko et ^l2004[2006 ). [^aarashvili et all 


2 OI 2 I deduced that there is no cut-off frequency for Alfven waves, only a damping 


mechanism due to the ion-neutral interaction effects, when one is taking into ac¬ 
count both the inertial and Hall term (where the inertial term is defined as the time 
variation of the relative velocity perturbations perpendicular to the unperturbed 
magne tic field between ions and neutrals). However, according to IVranies fc Kon^ 
(1201411 . this result might be due to the assumption of neglecting the exchange of 
momentum between ions and neutrals. However, one must be really careful with all 
these results since the ion-neutral collision frequency has been calculated in a very 
simplified manner and a proper treatment of the collision between ion and neutrals 
may provide different results (see below). In fact, even in semi-empirical models, 
using different methods to calculate the ion-neutral collision frequ ency provide very 


differ ent values for the damping and ambipolar diffusion terms (iDe Pontieu et al 


2001 1 . 


_ A mbipolar diffusion allows magnetic field t o diffuse into the atmosphere (ILeake fc Arber 

I2OO6I : lArber et al. 2007 : Leake fc Lintonl[201,'lh . Since neutrals move through the ex¬ 
panding magnetic field in the chromosphere, the magnetic field lines lift less matter. 

As a result the ambipolar diffusion inhibits the Rayleigh-Taylor (RT) instability. 

In contrast, a completely different configuration, such as the 2.5D simulations of 
prominences including partial ionization, shows an increase of small scale veloci¬ 
ties as a res ult of the non-linearity of the RT instab ility, in contrast to single fluid 
simulations ( Diaz et al iKhomenko et al. 1 I2H1 . Note that, the ambipolar dif¬ 

fusion can have effects of opposite nature: it can either inhibit the RT instability 
or increase the small scale structures of the RT instability dependi ng on the spatial 
distrib utio n of ambipolar diffusion, i.e., compare Arber et al\ i 200^ with lPfaz et al\ 

(|2014li and IKhomenko et an (l2014r ). Since the ambipolar diffusion is highly depen- 
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dent on the thermal properties of the plasma, the ion-neutral collision frequency 
needs to be calculated directly from the thermal properties of the plasma and the 
energy balance must be solved self-consistently. 

Regardless of the assumptions of the background atmosphere and the calculation 
of the ion-neutral collision frequency, it is clear that partial ionization effects play 
an important role in these processes. 

As mentioned before, ion-neutral interaction effects cannot change the topology 
of the magnetic field, but they can impact the reconnection rate. On one hand, 
ambipolar diffusion causes sharp structures to evolve around the reconnection re¬ 
gion, thus setting the scene for tearing reconnection, accelerating the reconnection 
rate or, depending on the configuration of the magnetic field and the diffusivi- 
ties, singularities that initiate reconnection. On the other hand, electron pressure 
and ohmic dissipation act against the formation of sharp currents, i.e., the re- 
connection rate depends on the ioniz ation rate and ion-neutral collision frequency 
( Brandenburg fc Zweibell 1994 . 1995 ). When the ionization decreases , the recon¬ 
nection rate decreases too ( Smith fc Sakai 2008 : Sakai &: Smith 2009ll . Note that 
the two fluid simulations from the latter two references used fixed ionization frac¬ 
tions and did not take into account the ionization/recombination rates. In addi¬ 
tion, they did not clarify if their parametric study is for a fixed density or a fixed 
number of ions since an increase of density can lead to a decrease of the recon¬ 
nection rate. In a reconnection process in a weakly ionized plasma when ions and 
neutrals are decoupled, it is found that an excess of ions c an accumulate in the 
reconn ectio n region, as shown wi th the simulations done by IVishniac fc Lazarian 
( I 999 II and Lazarian et al. 1 2004h . The recombination in the reconnection region 
removes ions which combined with the alfvenic flows can lead to fast reconnection. 
This has been tested with 2D two fluid simulations taking into account the ion- 
ization/recombination, op tically thin radiative transfer and collisional heating by 


Leake et all (l2012l . l2013a) . An important remark from the latter authors is that it 


is crucial to treat the ambipolar diffusion while at the same time solving the energy 
balance including the ionization/recombination since it can change drastically the 
impact of the ambipolar diffusion on the reconnection process. While the results 
described above indicate that ion-neutral effects are an important ingredient in a 
variety of processes in the solar atmosphere, the current simplified approaches to 
the energy balance significantly limit the applicability of these results in the solar 
atmosphere. 

All the studies mentioned in this section are focussed on specific driving pro¬ 
cesses. These studies are extremely useful to discern the physical mechanism in 
each specific problem. However, in order to place these processes in the context 
of the solar atmosphere, we need to combine this knowledge with self-consistent 
radiative MHD models that have the proper convective drivers (see Section [2|) and 
that inc lude the relevant phy si cal processes i n each layer of the atmosphe re, as 
done by Hansteen et al. ( 2007 1. AbbettI ( 2007 1. Martfnez-Svkora et all (l2C)08ll . and 


Hansteen et aL~(l2010l l. amongst others. The missing ingredient in these models 


has so far been ion-neutral interactions whic h should be calculated sel f-consistently 
with the thermal properties (such as done bv iLeake et aLl[^0I2L 12013^ for simplified 
magnetic reconnection simulations). The Hall effect in the photosp here has been 
implemented follo wing this approach by Cheung fc CameronI ( 201211 . The prelimi¬ 
nary work done in iMartfnez-Svkora et all (120121) added ambipolar diffusion to the 
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description of the outer solar atmosphere. We describe this type of simulations 
including ambipolar diffusion in the next sections. 


(c) Importance of the Ambipolar Diffusion in the Chromosphere 


In the previous section we described several aspects of how ion-neutral interac¬ 
tion impacts various physical processes in the chromosphere and prominences. Let 
us now discuss the range of values for the various “diffusion” terms of the expanded 
i nduct ion equation in the lower solar atmosphere. Figure 1 in iKhomenko &: Collados 
(|2012ll shows, u sing the semi-empirical VAL-C mod el and thin magnetic flux tube 
approximation ( SpruitllTQSlt Roberts fc WebbI 1978ll . the values of the various “dif¬ 
fusion” terms through the atmosphere. In short, ohmic diffusion dominates below 
the photosphere. In the photosphere, the Hall term is the largest and reaches 7 x 10® 
m^ s“®. In the chromosphere the largest, by roughly 2 orders of magnitude, is the 
ambipolar diffusion (10® m^ s“®). However, these values become extremely different 
if we consider a radiat i ve MHD simulation such as th e one used by, for instance, 
Gudiksen fc Nordhmdl (12005 a b Hansteen et al ( 2010ll (see below). 


In order to calculate the ambipolar diffusion terms, the ion-neutral collision fre¬ 
quency must be obtained from the ionizati on state and thermal properties of the 


plasma. For this there are several metho ds (lOsterbrockI 1961 


_ _ _ ----- - - _ . _ _ - Ivon Steiger fc Geisa 

19891 : Fontenla et al\ 199, ll among- others) . OsterbrockI ( 196lh uses el astic scattering 


cross s ections . Ivon Steiger fc GeissI (1198911 considers hard spheres . andiFonteiila et al 
(1993) uses different sets of elastic collision cross sections than OsterbrockI ()l96i l 


to determine the ion-neutral collision frequency. The final results of these methods 
differ in the ion- neutral collision frequenc y by an order of magnitude for semi- 
empiric al models ( De.Jk)ntieuel^ ( 200lll among others). In a more realistic ap¬ 
proach, lMartfnez^Sytor^.eLa^ ( 2 OI 2 II show that the various methods used to calcu¬ 
late the ion-neutral collision frequency differ most in chromospheric regions where 
the ambipolar diffusion is l argest. A more recent a nalysis of the ion-neutral col¬ 
lision frequency is done bv IVranies fc Krstic ( 2013 1 where they include accurate 
cross sections which vary as a function of temperature, due to quantum effects, 
and between descriptions of elastic scattering, momentum transfer and viscosity. 
As mentioned by these authors, this work can be improved including the effects of 
inelastic collisions. 

Figure [1] show s the ohmic diffusion, Hall term and ambipolar diffusion where 


OsterbrockI (Il96ll l was used to calculate the ion-neutral collision frequency. The 


simulation is 2.5D and ranges from the upper convection zone to the lower corona 
i ncluding radiative transfer with scattering fr om the photosphere to the corona 

( Havek et al.ll2010l : Carlsson fc Leenaartsll2012 1. thermal conduction a long the mag- _ 

netic field lines and partial ionization effects using the Bifrost code ('see lMartfnez-Svkora et al. 
I2OI2I and references cited within). The initial magnetic field is unipolar with a mean 
unsigned field strength of 5 G in the corona. The convective m otions build enough 
magnetic field stress to self-consistently maintain the hot corona dGalsgaard fc Nordhmdl 


1996 : ^diksen fc Nordhmd200412005 1)1: iHansteen et Qi.ll2010l : lMartmez-Svkora et al\ 


2011all . Note that the magnetic field is vertical in the corona (left panel), and despite 
the fact that our simulations do not include imposed flux emergence, the dynam¬ 
ics in the simulation and the convective motions accumulate horizontal magnetic 
field in the sub-adiabatic photosphere. As mentioned above, this horizontal field 
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will in principle be c onstrained to the sub-adiabatic photosphere (Acheson 1979t 


Archontis et al.l[2004 1. In addition, magnetic field footpoints will move around. Note 


that these simulations naturally reproduce several of the energy generating drivers 
discussed in Section [2] Moreover, the energy flux through the atmosphere will be 
modified by the physical processes considered in the model as it propagates up¬ 
wards towards the corona as mentioned above Isee lMartinez-Sykora et a7l2012 and 
compare with Section [3]). 


t; [mVs] 



0 5 10 15 0 5 10 15 0 5 10 '5 

« [Mm] * [Mm] . [Mm] 


Figure 1. Maps of ohmic diffusion (left). Hall term (middle) and ambipolar diffusion (right) 
from a snapshot of a 2D radiative-MHD simnlation that includes partial ionization effects. 
The color scheme is on a logarithmic scale. The magnetic field is drawn with white lines 
in the left panel. 


The ambipolar diffusion is the largest term in the induction equation in the 
chromosphere by almost ten orders of magnitude compared to the ohmic dissipa¬ 
tion and something b etween two and six orders of magnitude larger than the Hall 
term (see Figurellland lMartmez-Svkora et alll2012l the latter compares these terms 
with the hyper-diffusion intrinsic in the code too.). In addition, the various terms 
show large spatial variability. For instance, the ambipolar diffusion changes at very 
small scales (within a few hundred km) by up to seven orders of magnitude, e.g., at 
height z = 1 Mm. This large variability is due to the large changes in the ion-neutral 
collision frequency and in the neutral density when the ionization state and thermal 
properties of the fluid are taken into account. Note that this model is self-consistent 
and the ionization state and thermal properties are set by the various drivers, such 
as shocks, and convective motions of the footpoints that build magnetic field stress 
(Section [2|) and the physical processes implemented in the simulation: radiation, 
thermal conduction along the magnetic field lines, scattering, partial ionization ef¬ 
fects in the equation of state and ion-neutral interaction effects. One important 
process that has yet to be taken into account in these models which will change 
the ohmic diffusion, Hall term, and ambipolar diffusion is the effect of time de¬ 
pendent non-equilib rium ionization of hydrogen and helium ( Leenaarts et al ] l2nn7l: 


Golding et aZ.ll2014D which are know to be important in the upper chromosphere. In 


our models we assume time dependent equilibrium ionization. These effects must 
be included in order to draw firm conclusions on the importance of ion-neutral 
interactions in this part of the solar atmosphere. 

From Figure[I] it is clear that any magnetic energy flux (driven by photospheric 
or chromospheric processes) that is transferred through the chromosphere may be 
strongly influenced and/or diffused by the ambipolar diffusion. We can therefore 
already conclude that the chromosphere may diffuse more magnetic field stress 
than any other atmospheric layer. Note that the ohmic diffusion is greater in the 
photosphere and chromosphere than in the corona. 
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Table 1. Description of the simulations 

(The left column lists the names of the various 2D simulations, and the right column 
gives a short description of each simulation.) 


Name 

Description 

NGOL 

without ion-neutral interaction effects 

GOL-OS 

with ion-neutral interaction effects using Osterbrock llQB!! for 

GOL-F 

with ion-neutral interaction effects using Fontenla et al. llQQSl for Uir, 


(d) Impact on the Solar Atmosphere 

Here we will not describe in detail how ambipolar diffusion affects the various 
processes such as flux emergence and wave propagation in 2D radiative MHD sim¬ 
ulations. This is outside of the scope of this review and will be described in a series 
of follow up papers. However, we will show the impact of ambipolar diffusion on 
the thermal properties of the solar atmosphere. For this we compare different 2D 
radiative-MHD simulations listed in table [TJ 

The different models, i.e., NGOL, GOL-OS and GOL-F, show important dif¬ 
ferences. This means that even with the large artificial diffusion intrinsic in these 
codes, the resulting models can discern (at least) some of the effects produced by 
the partial ionization interaction effects. A first look at the 2D temperature maps in 
Figure [5] reveals at least two thermal properties that differ between the various sim¬ 
ulations: 1) the cold chromospheric expanding bubbles have higher temperatures in 
the GOL simulations, and they do not rise as high as in simulations without partial 
ionization effects. 2) The transition region shows a smaller temperature gradient 
(i.e., is spread out over a larger height range) in the GOL-OS simulation than in 
the other two simulations. 


Temperolure [k] time=8200s 


2-OE+D3 3.6E+03 6.3E+03 1.1E+04 2-0E+04 



Figure 2. Temperature maps for the NGOL simulation (left), GOL-OS (middle) and 
GOL-F (right) reveal differences in the thermal properties. The temperature is shown 
in logarithmic scale. 


Gold chromospheric bubbles are formed by expanding shock fronts that pass 
through the chromo sphere. These bubbles contain the lowest temperatures in the 
atmospheric model (jLeenaarts et al.N20111l . Statistically, the simulations with par¬ 
tial ionization effects have hotter and denser cold chromospheric expanding bubbles 
than the NGOL simulation. This is shown in the Joint Probability Distribution 
Functions (JPDF) of the density and temperature integrated over 30 minutes in 
Figure |3] for simulations NGOL (left panel), GOL-OS (middle panel), and GOL- 
F (right panel). Everything below 4000 K (logT < 3.6) corresponds to the cold 
chromospheric bubbles. The minimum temperature achieved is lower in the NGOL 
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simulation than in the GOL simulations (enhanced with white circle E in the middle 
panel). In fact, the coolest temperatures for the NGOL simulation are controlled 
by an ad-hoc heating term in the Bifrost code which is introduced to avoid reaching 
temperatures lower than ~ 1600 K since below this threshold the tabulated equa¬ 
tion of state is no longer accurate ( Garlsson fc Leenaartsl[2012 1. While this ad-hoc 
heating is crucial for the NGOL simulation, it is rarely necessary in both simula¬ 
tions with ion-neutral interaction effects since the cold chromospheric bubbles are 
not cold enough to reach the temperature at which the ad-hoc heating starts. The 
cold chromospheric bubbles in the simulations with partial ioniz ation effects are 
a few hundred of degrees hotter than in the NGOL simulation. iLeenaarts et a,I 


(|201l[l found no mechanism to prevent extremely low temperatures from occurring 
in expanding bubbles within the standard, single-fluid MHD formulation. The joule 
heating contribution coming from ambipolar diffusion provides such a mechanism. 
Observations can help dete rmine the amount and temperature of s uch cold plasma 
in the solar chromosphere ( Avres fc RabinI 1996; Penn et al. 2011 1. The details of 
the physical processes in the chromospheric cold bubbles will be described in a 
follow up paper. 


Log(Histogram) 


'V.8E-12 4.8E-09 4.aE-06 4,eE-03 4.8E+00 



-16 -14 -12 -10 -e -16 -14 -12 -10 -8 -16 -14 -12 -10 -8 

Log(Derslty) [g/cm^] LQg(Density) [g/cm^] Log(Density) [g/cm^] 


Figure 3. Joint Probability Density Function (JPDF) of temperature (vertical axis) and 
density (horizontal axis) integrated over 30 min for simulations NGOL (left), and GOL-OS 
(middle) and GOL-F (right). The white contours correspond to the temperature and 
density regime of the simulation NGOL in dotted line, GOL-OS in solid line and GOL-F 
in dashed line in order to simplify the comparison between the simulations. The circles in 
the middle panel enhance the differences between the various simulations. 


The upper chromosphere also shows considerable differences between the three 
simulations. One can appreciate that in the GOL simulations there are a greater 
number (compared to NGOL) of grid points concentrated around the temperature 
where hydrogen ionizes in the simulations (white circle D in the middle panel of 
Figure [3]). Moreover, the number of grid points around the temperature where he¬ 
lium ionizes {T ~ 10"'’ K) is greater for GOL-OS compared to the other simulations 
(circle C). In fact, GOL-OS shows several grid points with rather dense material 
in the transition region (circle A) while the transition region is the least dense in 
NGOL (circle B). GOL-F has the narrowest range of densities in the transition re¬ 
gion. These discrepancies in the GOL simulations are due to the differences in the 
distribution of ambipolar diffusion as a result of the different methods to calculate 
the ion-neutral collision frequency and cross section values, hydrogen and helium 
time dependent ionization may spread the amount of plasma over a larger range 
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of temperatures than for the case of ionization equilibrium. This latter case leads 
to the temperature of the plasma being localized within the values where hydrogen 
(circle D) and helium are ionized (circle C). 

There are several physical processes that lead to these differences in the upper 
chromosphere. These processes will be described in detail in a follow up paper. In 
short, the expanding cold bubbles, which have large ambipolar diffusion, are heated 
and, it is through these bubbles that the photospheric horizontal magnetic field is 
diffused into the chromosphere. From there, the field expands into the corona push¬ 
ing chromospheric material to higher layers. A large percentage of the magnetic 
stresses due to this expansion of the magnetic field, waves and footpoint motions 
are dissipated in the cold bubble. Moreover, for GOL-OS, the diffused field lines 
coming from the photosphere interact and dissipate magnetic energy in the upper 
chromosphere into thermal energy driving temperatures up to ^ 10'^ K. The vari¬ 
ous regions of the atmosphere will impact how mechanical and magnetic energy is 
transported to greater heights in different manners, since different processes domi¬ 
nate in the various regions. The partial ionization interaction effects also produce a 
slightly more dynamic upper-chromosphere, transition region and lower corona as 
a result of the expansion of the diffused magnetic field and the interaction of this 
expanding magnetic field lines with the ambient magnetic field. These simulations 
need to be r ecalculated taking into account a better ion-neutral collision frequency 
calculation ( Vranies fc Krstic 201, “ih . Including the effects of partial ionization is 
important to many aspects of the physics of the chromosphere and corona, in par¬ 
ticular those related to the evolution of the magnetic field and the energization of 
the outer atmosphere’s plasma. 


5. Summary 

Including the effects of partial ionization is important to many aspects of the physics 
of the chromosphere and corona, in particular those related to the evolution of the 
magnetic field and the energization of the outer atmospheres plasma, as the inter¬ 
action between ionized particles and neutral particles has important consequences 
on the magneto-thermodynamics. In particular, our simplified 2D radiative MHD 
models reveal that the ambipolar diffusion dissipates magnetic energy and increases 
the minimum temperature in the chromosphere and the upper chromosphere may 
get heated and expanded over a greater range of heights. 
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